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We analyse the influence of a surface viscosity on the orientational 
dynamics of a nematic liquid crystal subjected to an oscillatory Cou- 
ette flow. Approximate analytical solutions of nematohydrodynamic 
equations for small flow amplitudes are calculated and compared with 
the results of full numerical simulations. The range of flow frequen- 
cies where the surface viscosity has strong influence on the optical 
response is determined. 
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INTRODUCTION 

In a nematic liquid crystal (NLC) complex flow behavior arises from 
the strong coupling between the velocity and orientational degrees 
of freedom (director n). Orientational dynamics of NLC under oscil- 
latory flow has generally been studied in the case of strong surface 
anchoring (fixed orientation at the substrates). The influence of weak 
surface anchoring and surface viscosity on the dynamics of director 
reorientations under an applied electric field has been discussed in 
0, [|. The surface-dominated orientational dynamics in NLC em- 
bedded in a solid porous matrix has been studied by means of the 
dynamic scattering technique, from which the temperature behavior 
of surface viscosity has been deduced ||. Recently, a flow-induced 
orientational transition in nematics at the substrate with weak planar 
anchoring was found and an estimate of surface viscosity has been 
given g. 

In this paper the orientational behavior of the NLC layer with 
weak surface anchoring subjected to the rectilinear oscillatory Cou- 
ette flow is studied theoretically. We analyse the case of small flow 
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amplitudes where the director motion is within the flow plane. The 
influence of anchoring strength and surface viscosity on the orienta- 
tional dynamics are investigated and the appropriate experimental 
conditions to measure the surface friction are proposed. 



BASIC EQUATIONS 

The nematic layer of thickness d is confined between two identical 
substrates which provide weak surface anchoring at z = ±d/2 in a 
Cartesian coordinate system. The oscillating flow is in the x direc- 
tion and the director confined to x — z plane. In this situation the 
director and the velocity are only functions of the distance z from 
the boundaries and time t 

n x = cos 0(2, t), n y = 0, n z = sxn9(z,t), 

v x = v x (z,t), v y = 0, v x = 0, (1) 

where 9(z,t) is the angle with respect to the x axis. With the di- 
mensionless variables z — z/d, t — tu, v x = v x /(ud) the equations 
governing the alignment and the flow can be rewritten as [g], §] 



d >t - K(9)v X;Z = e 



P(0)9, zz + \p' 



(2) 



Sv Xit = d z {-{\ - \)K(d)9, t + Q(9)v x , z }, (3) 

where the tildes have been omitted and 

A cos 2 9 - sin 2 9 a 3 2 . 2 K 33 

K{9) = , A = — , P[9) = cos 9 + K sin 9, K = , 

1 - A a 2 An 

2(— oii)Q{9) = «4 + (ocq — CK2) sin 2 + (a 3 + a6 + 2ai sin 2 6 1 ) cos 2 9, 

^ lid 2 pd 2 

-,r<i = -r^,7i = a 3 - a 2 ,o = r v u,r v 



T d UJ An -«2 

Here the «j are the viscosity coefficients, iC^ are elastic constants, p 
is the density of NLC and the notation = df/di, f'(g) = df/dg 
has been used throughout. Boundary conditions for the velocity for 
the oscillatory Couette flow are 

v x (z = -1/2) = 0, v x (z = +1/2) = acos(t), (4) 

where a = A x /d and A x is the upper plate displacement amplitude. 
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The weak anchoring is described mathematically in terms of sur- 
face energy per unit area F s = (W/2)f s {6 — 6> ), where W is the 
anchoring strength and function f s {9 — 9q) has a minimum at 9 = 6q. 
A simple phenomenological expression for the surface energy was in- 
troduced by Rapini and Papoular where f s (9 — 9 ) = sin 2 (6> — 9 ) |7|]. 
The boundary conditions for the director can be obtained from the 
surface torques balance equation |2j, || 

Wd df a ood 39 n 

— on z = — 1/2 and + on z = 1/2. 

Here rj = 7i/ 71 is so-called surface viscosity which characterises the 
dissipation at the substrate; 71 is the bulk orientational viscosity and 
l ll is a characteristic interfacial length for surface viscosity. 

In order to obtain an approximate analytical solution of Eqs.(H), 
(D) we consider the case of small flow amplitudes a < 1 which corre- 
sponds to the small distortions of the director profile 

9 = 9 + 9, v x = v x0 + U, |0|<1, \U\<^1, (6) 

where #0 = const, v x o = is the solution of (||), @ in the absence 
of oscillatory flow (a = 0); the value #0 is defined by the minimum of 
surface energy given for small amplitude of the director oscillations 
as / s = 9 2 . In the low-frequency range to be considered here one has 
5 <sC 1 (p pa 10 3 kg/m 3 , \a 2 \ ~ 1CU 1 N-s/m 2 and d « m gives 
5 < 1 for frequencies / < 1 kHz) and the inertia term [left-hand side 
of Eq.(|3|)] is neglected. Then one obtains from equations @, (Q) 

9,t - K U, Z = eP 9, zz , 

(1 - X)K 9 jtz - Q U ZZ = 0, (7) 

where K = K(9q), Po = P(9o), Qo = Q(9 ), with the boundary 
conditions 

9, z — E9 — G9 )t = |z=_i/2 , 9 >z + E9 + G9 )t = \ z =i/2, 
U(z = -1/2) = 0, U(z = 1/2) = acos(t). (8) 

Here E = Wd/(P ^n), G = ud^J (PqK^). Since Eqs.© are 
linear, the periodic boundary conditions (H) (periodic forcing) will 
lead to time-periodic solutions for the director 

9{z,t) = Ti(z) cos(t) + T 2 (z) sin(t) (9) 
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with 



T^z) = aK ( 



o 



C x F x {z) - C 2 F 2 {z) 
C\ + CI 



, T 2 (z) = -aK 



C 1 F 2 (z) + C 2 F 1 (z) 

ci + Cf 



and for the velocity 



U(z,t) 




(10) 




2Q k C\ + C\ 



where the functions F^z) and Q are given in Appendix A. 
RESULTS AND DISCUSSION 

For weak homeotropic anchoring (6 = it/2, W = 10~ 6 J/m 2 ) the 
profiles of the director deviation Ti(z), T 2 {z) from the homeotropic 
position as well as the velocity components Ui(z), U 2 (z) are pre- 
sented in Fig. 1 for different values of the surface viscosity i] = 71 / 7l . 
The material parameters of MBBA (see Appendix B) have been used 
throughout. The curves for strong homeotropic anchoring [fixed ori- 
entation at the substrates, 9{z = ±1/2) = 0] as well as for zero 
surface viscosity (Z 7l = 0) are also shown for comparison. For the 
velocity component U\ we subtract the linear part a(z + 1/2) which 
corresponds to the velocity profile in case of isotropic liquid. One 
can see that the value of surface viscosity strongly influence on the 
director and velocity distribution under the oscillatory Couette flow. 
Note, that raising of the surface viscosity has a similar effect as the 
increasing of the surface anchoring strength W. 

The amplitude of the director oscillations is proportional to Kq = 
K(6 ). For MBBA one has K(n/2)/K(0) w 10 2 and the amplitude 
of the director oscillations for weak planar anchoring is much smaller 
than for the homeotropic orientation at the same flow amplitude a. 
Therefore, in case of planar orientation the influence of the surface 
viscosity on the orientational dynamic of NLC is much smaller than 
for weak homeotropic anchoring. 

One of the widely used technique for studying the orientational 
behavior of NLC is the measurements of the transmitted through 
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Figure 1: Director and velocity profiles for oscillatory Couette flow. 
a = 0.2, / = 5 Hz, d = 10 fim, W = lO" 6 J/m 2 , Z 7l [m]: (— ); 10~ 7 
( — ); 10~ 6 (•••), strong anchoring (-•-). 



the NLC cell polarized light. In the geometry where the polars are 
crossed and the x axis is at 45° one has for the light intensity 



1/2 
-1/2 



n n e 

= - n. 



2 cos 2 + n 2 sin 2 



dz , (11) 



where n G , n e are the ordinary and extraordinary refractive indices 
and A is the wavelength of light. Time-periodic director oscillations 
lead to a variation of the optical response I = I(t). Using the di- 
rector solution (H), the transmitted light intensity can easily be cal- 
culated. We found that depending on the value of surface viscosity, 
the maximum of the intensity changes and its position is shifted in 
time with respect to the moment of zero displacement of upper plate 
(t — 7r/2). In order to find the range of flow frequencies where the 
surface viscosity has a strong influence on the optical response, the 
dependence of the maximum of the transmitted light intensity I ma x 
on the flow frequency was calculated for different values of surface 
viscosity (Fig. 2). We present the results for two typical values of 
anchoring strength W within the range of recently found experimen- 
tal ones . The optical response for the case of strong anchoring is 
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Figure 2: Maximum of the transmitted light intensity versus the 
frequency of oscillatory Couette flow at a = 0.2. (a) - d — 10 /zm, 
W = 10~ 6 J/m 2 ; (b) - d = 5 /mi, W = IO" 5 J/m 2 . Z 7l [m]: (— ); 
10~ 7 ( — ); 10~ 6 (•-■), strong anchoring (-•-). 



shown for comparison. It is clearly seen that there exists some flow 
frequency range [~ 0.1 4- 10 Hz in Fig. 2(a) and ~ 1 4- 100 Hz in 
Fig. 2(b)] where the optical response is strongly depending on the 
value of surface viscosity. Below this frequency range the transmit- 
ted intensity is small and the determination of the surface viscosity 
is complicated, whereas at higher frequencies the boundary layers of 
the order of J 'K u / '(2uji) at z — ±1/2 |J become thinner and the 
influence of the surface viscosity on the director dynamics is not so 
apparent. 

In order to verify the approximate solutions, direct numerical sim- 
ulations of system (fj), @ with boundary conditions (H), @ have 
been performed. It was found that for the frequency range u t^ 1 
(i5< 1) one can safely use the analytical small- amplitude solutions 
up to the flow amplitudes a pa 0.25 for the homeotropic anchoring 
(the difference between analytical expressions and numerical solution 
does not exceed 0.5%). 
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In conclusion, the influence of the anchoring strength and surface 
viscosity on the director dynamics in the bulk of NLC and at the sub- 
strates is investigated for oscillatory Couette flow. Approximate an- 
alytical solutions for small-amplitude oscillatory flow in case of weak 
surface anchoring are obtained. The analysis shows that there exists 
a certain flow frequency range where the optical response strongly 
depends on the value of surface viscosity. This effect can be used 
for a more precise experimental determination of the surface viscos- 
ity and investigations of surface-dominated orientational dynamics in 
oscillatory flow of NLC. 



Appendix A: 



Fi(z) = cosh(A;2;) cos(kz) — cosh(/c/2) cos(fc/2) + 

+ ^ 2 fc ^ 2 {cosh(A:/2) sin(fc/2) - sinh(A;/2) cos(A;/2)} - 
hG 

- ^ 2 + G , 2 {cosh(A:/2) sin(fc/2) + sinh(A;/2) cos(A;/2)}, 

F 2 (z) = sinh(kz) sm(kz) — sinh(/c/2) sin(/c/2) — 
kE 

r {cosh(A;/2)sin(A;/2) + sinh(A;/2)cos(A;/2)} - 



E 2 + G 2 
kG 



{cosh(A;/2) sin(fc/2) - sinh(A;/2) cos(A;/2)}, 



E 2 + G 2 

F 3 (z) = sinh(A;2;) cos(kz) — cosh(kz) sm(kz) — 

-2z(smh(k/2) cos(A;/2) - cosh(A;/2) sin(fc/2)) , 

Fi(z) = sinh(A;2;) cos(kz) + cosh(fcz) sm.(kz) — 

-2z{smh{k/2) cos(A;/2) + cosh(A;/2) sin(fc/2)) , 



k 



;i - mi 



2eQ P 
d = sinh(A;/2) sin(fc/2) - 
(1 - \)K 2 ' 



Qok 
kE 



cosh(A;/2) sin(fc/2) - sinh(A;/2) cos(A;/2) 



{cosh(A;/2) sin(Jfe/2) + sinh(A;/2) COS [K 



+ 



E 2 + G 2 
kG 

, ^ n {cosh(A:/2) sin(fc/2) - sinh(A;/2) cos(A;/2)}, 



E 2 + G 2 
C 2 = cosh(A;/2) cos(A;/2) 



s 
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(I _ X\K 2 

- cosh(fc/2) sin(fc/2) + sinh(fc/2) cos(A;/2) 

QqK 

--= T=r{cosh(fc/2) sin(fc/2) - sinh(fc/2) cos(£;/2)} + 

E 2 + G 2 J 

+ E2 + G2 {cosh(fc/2) sin(A;/2) + sinh(A;/2) cos(£;/2)} . 
Appendix B: 

The numerical computations are carried out for the following MBBA 
material parameters at 25 °C || [10]: viscosity coefficients in units 



of 10" 3 N-s/m 2 : oti = -18.1, a 2 = -110.4, a 3 = -1.1, a A = 82.6, 
as = 77.9, «6 = —33.6; elasticity coefficients in units of 10~ 12 N: 
K n = 6.66, K 2 2 = 4.2, K 33 = 8.61; mass density p = 10 3 kg/m 3 and 
refractive indices for wavelength of light A = 670 nm: n = 1.542, 



n e = 1.7435 UTT 
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